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ABSTRACT 


A flight loading profile based solely on accelerometer 
data is developed ini cl AcootaES for all fatigue damaging 
loads experienced by high performance aircraft. A random 
sequence of the flight loads is generated in a computer 
program thet applies the ESDU Cumulative Damage Hypothesis 
fOr Gamage prediction of fleet aircraft. Data from actual 
accelerometer readings are used in evaluating whether 
reasonable results can be obtained from the program. The 
feasibility of implementing the program is discussed along 
with recommendations improving the suitability of the 


program for more accurate predictions. 
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Pe Se NCReDUC EON 


Mies designe phaser of new adrcraft stnvolves#the use of 
an estimated flight loading profile in order to make a ser- 
vice fatigue life estimate of the aircraft. A full-scale 
fatigue test is normally performed on the aircraft in the 
early stages of production to verify the usable life predic- 
tion; however it 18 inevitable that the actual loads the 
aircraft experiences throughout its service life will be 
different than the design spectrum. The fatigue monitoring 
program presently used for naval aircraft attempts to record, 
without regard to sequence, the actual loads an aircraft 
experiences by means of a four level exceedance accelerometer. 
Anabyseismoterilight data from oscillograph recordings indicates 
that each type of aircraft has a culicineatshabile flight load- 
Ineeeeoetmleeond i¢ has bDeem determined that the order in 
which these loads occur on a particular type of aircraft 
hassamsipnificant effect on the damage prediciton. This 
thesis proposes a procedure which estimates the actual flight 
profile from accelerometer data and statistically formulates 
a probable order in which these loads may have occurred. A 
fatigue cumulative damage theory is then incorporated wnich 
accounts for the sequential effects of loads in fatigue damage. 
The method is adapted to the computer which is necessary for 
monitoring the fatigue damage of a Jarge fleet of aircraft 
and provides a damage calculation for each aircraft in 


meron mMately 2.5 seconds. 
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II. DAMAGE PREDICTION 


Structural components which experience repeated loadings 
are susceptible to fatigue damage for loading amplitudes 
which occur well below the fracture strength of the material. 
Mneses LOdcdine@ yecondi tions oecur frequently on aircraft, since 
woe requirement of maintaining a particular flight pattern 
requires many loading changes. As a result, the design of 
aircraft structural components must account for the damaging 
effects of these repeated loadings, which necessitates an 
accurate means of predicting the fatigue damage accumulated 
Gue to these loading conditions. 

Fatigue life will be understood to be the life until a 
visible crack is present, or the life until complete failure 
of a small component where the crack propagation phase 
WOlmG@mbDemshort in compariison to the crack initiavion period. 
Damage will refer to the percent of fatigue life expended 
due to a prescribed loading condition. Methods for predic- 
ting damage invariably rely upon some application of the 
time-tested Palmgrem-Miner Rule, better known simply as 


Miner's Rule. 


Aj MINER S&S RULE 

Both Palmgren and Miner are credited with introducing 
what has become the most well known Incremental damage 
reory (Ref. 1). If a component has experienced n, load 


ioud amplitude, the 


° 
« 


cycles with the same mean load and 





portion of the fatigue life expended is n, /N where N is the 
life to failure in a constant-amplitude test with the same 
mean and amplitude. Failure is assumed to occur if the sum 
of the consumed life portions equal 100 percent. MThis 


ampiies tuhatethe condition for failure? ils: 
mn, /N = (ip) 


This failure criterion was extended to a component 
experiencing load cycles of various amplitudes and means 


PpUCtM Chale EMe COnd1tion for failure is: 
Z(n,/N, ) = il C25 


Specimen testing of different materials under varying 
mean and mrt ude loads provide the N. to use in Miner's 
Rule. The family of curves developed from specimen testing 
for a particular stress concentration factor are referred 
to as S-N curves. 

itm Advantages 

Due to the simplicity of application, Miner's rule 
has enjoyed widespread popularity as a means of predicting 
Solace Sin Jo brorouGae 1% 1Ss Darbicularly suited to.a large 
fleet of aircraft where loading data obtained from acccler- 
ometer readings can be readily stored in computers, which 
glows for fully automatic processing of data and evaluation 


meedamace prediction for individual aircraft. Mtner'’s rule 


eG 





provides a quick "rule-of-thumb" indication of fatigue life 
expended since as the number of cycles of loading increases 
a greater portion of fatigue life is expended. 

Om LAmMLcations 

Due to the complexity of aircraft loading cycles, 

Miner's rule has proven to be a rather conservative estimate 
of fatigue damage. For instance Lambert [Ref. 2] applied 
a tuyelcal aircraft block program to a rivetted lapjoint of 
2024-T3 clad sheet aluminum. The ratio cf actual life to 
Miner's estimated life obtained was 2.64. Through an exten- 
sive survey of the area of fatigue damage Schijve [Ref. 1] 
concludes that a Miner's Rule approximation provides a con- 
siderably @maccuratey estimation of the fatigue life expended 
bye oli eter CONOM CcollyeaT 1S nol isound to retire an 
airerait trom Service when actually only 37 percent of its 


fatigue life has been expended. 


B. SEQUENCE EFFECTS 

Since loading amplitudes vary randomly on an aircraft 
Ue tOlUuum US lic, stuemas been determined that vhe jrder 
in which these loads are applied must be significant in 
the fatigue prediction. S-N data applied to Miner's Rule 
are based on data points obtained by specimen failure at a 
CGmeilame arpiatude for a’ specific mean stress; therefore, 
mois higily Unlilely that damage accumulation in aircraft 
components, which experience a wide range of amplitudes from 


eeele, co cycle, would accurately adhere to a Miner's Rul 


al 





approximation. The order of loads is not taken into account 
by Miner's Rule. Schijve [Ref. 1] attributes the difference 
in damage accumulation between variable fatigue loading and 
constant amplitude loading to interaction effects. 
1. Residual Stresses 
If a load amplitude is of sufficient magnitude to 
cause plastic deformation in the region of a crack tip, re- 
Sidual stresses will remain in this region unless the defor- 
mation is fully reversed. These residual stresses have to 
be added to the Lueees induced by the applied loads and 
as a consequence can significantly affect the fatigue damage 
accumulation. 
2. Cnack Closure Theory 
If a cracked sheet is loaded in compression the 
crack will be closed. Elber [Ref. 3] observed that crack 
closure may occur while the sheet is still loaded in tension. 
According to Elber,plastic elongation will occur in the plas- 
TlGma@ene Of bie growing crack. Thils' plastic deformavion 
will remain present in the wake or the crack and will cause 
closure before complete unloading of the specimen. 
Peep cavlon 
Highly maneuverable aircraft experience a number of 
high g loadings followed by a large number of smaller positive 
Paola danrSwee rior to completion of crack initiation thege 
high g loadings produce a residual stress at points of high 
stress concentration which effectively lower the mean values 


of the subsequent larger number of smaller positive g loads 





and consequently reduce the actual fatigue damage. Once 
crack initiation is complete, these high g loadings produce 
residual stress in the crack tips, which tends to retard 
the propagation of the crack and consequently reduce the 
actual fatigue damage. If an extreme negative loading 

were experienced, the residual compresSive stress induced 
would have the opposite effect and tend to cancel the favor- 
able effects from the extreme positive loading. Figures 
(la and 1b) demonstrate the results Schijve [Ref. 1] obtained 
from subjecting specimens of a known initial crack length 
to the different loading spectra. In Figure 1b note the 
significant number of additional cycles (point c) needed to 
obtain the same crack length after a few high positive 
loadings. The compresSive residual stresses induced by 
high loadings seems tooffer a plausible reason for the 

ia oie eter ative estimates of aircraft fatigue damage 
from a Miner's Bite approximation. It should be noted 
Chdiarone saomatingeeround—air—pround cycle tends to relieve 
the positive effects of extremely high maneuver load 
residual stresses and should be account for in fatigue 
damage prediction. — 

In view of the effect that load sequence has upon 
nate vempredi1cllon., 1t appears that not only are the loadings 
C=perlcnced Dy the aircraft imoortant in damage prediction, 
Mitcrakeo thie Order in which they are imposed. A later sec- 


tion will discuss the Engineering Science Data Unit (ESDU) 





Cumulative Damage Hypothesis, which takes into account 
sequence effects. The next section will address the prob- 
lem of estimating the actual load spectrum experienced by 


a given aircraft. 


A 





Tit. DATA ACQUISITION 


Pie transpOGihearltrcract boae Bust LOading spectrum is 
considered to be the most significant loading spectrum that 
contributes to fatigue damage. These loadings appear to be 
Gaussian in nature and occur symmetrically aoe the one-g 
level. Consequently, statistical methods such as the power 
spectral density are being employed to determine the spectrun 
of loads that these aircraft have experienced. However in 
high performance military aircraft, maneuver loads are con- 
Sidered the damaging loads for fatigue damage, and the gust 
loads are either disregarded or considered as a negligible 
centributor to fatigue damage. Leybold's attempt [Ref. 14] 
to describe the maneuver load spectrum from known probability 
Gis vempupton PuncttoOns met with some success’ but was re- 
Strlenco tO a particular type of aircraft. Data [Ref. 11] 
have shown significant differences in maneuver loads within 
aircraft of the same type as well as between aircraft of 
different types [Ref. 8]. The main reason for this differ- 
ence is the wide range of missions performed by these air- 
eraft and yet knowledge of the loads experienced by individual 
aireraft is necessary for accurate damage predictions. 

BiressCS ap tavieue critical points and their number 
of occurrences are required for fatigue damage prediction. 
Stress ae POlnu canioy ve recorded directly and conse- 


Meme vaiaranelLers tiatv relate to ctress mist be utilised. 


i 
ee 








The most widely used parameter for calculating stress 
is strain. The theory of elasticity has mathematically 
developed the relationship of stress to strain. Unfortunately 
PrewiOus, ALLEMDLS LO record strain history irom strain gauges 
have been unsuccessful due to the unreliability of gauge 
readings caused by zero drift and the difficulty in processing 
and reducing strain data for a large fleet of aircraft. 
Serateh gauges have also been recommended and are capable of 
PecCoradine many NOUrs O] Strain On @ tany disk but have the 
same logistics and reliability problems as strain gauges. 
Despite previous setbacks,strain recording devices are still 
considered to be the best means of obtaining stress. The 
National Aerospace Laboratory, Netherlands has developed 
a prototype system for in-flight strain recording with on- 
board data reduction Lice laeehius system themoulous som 
an ordinary resistance strain gauge bridge is amplified and 
digmerZzcamme nls difital record iS next analyzed for peak 
values using a "threshold" criterion. The values of succes~ 
sive minima and maxima are written on an on-board magnetic 
tape recorder. The prototype is currently installed in an 
operational F-10146 aircraft for long-term operational testing. 
Tries J res ullts are highly Dicom Same. 

Another method for determining stress,and the most widely 
Ueeemparant@per 10r calculating stress On aircraft,is acceler- 
ation. Depending on the model of the accelerometer, discrete - 
Levels of fs exceeded are digitally recorded. stress analy- 


PImommectioeiapeeue critical points provides a relationship 
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between acceleration and stress. This method is presently 


used by the U.S. Navy. 


A. COUNTING METHODS 

Actual load-time histories will consist of a number of 
load excursions with an irregular pattern and in an irregular 
sequenee. = The analysis of load-time histories has to be 
such that the amount of damage caused by these load excursions 
is quantitatively reflected in the final results. Reference 
5 describes ten different counting procedures and their 
somluty toerecord sitgenificant fatigue loads on aircraft. 

The methods are subdivided into three distinct categories 
depending on the parameters of the loading recorded. Basi- 
cally the counting procedures record either peaks (maximum 
or minimum values in Figure 2), level exceedances (for exam- 
ple levels one or two in Figure 2), or ranges (for example 
the distance between level one and minus one). The pre- 
ferred method of counting is the range-pair-range method 
where all ranges of Significance are counted. The drawback 
to this method is the extensive data reduction required. 

The restricted level~crossing count method is also con- 
sidered an effective means of recording fatigue damage. 
Accelerometers used in fatigue analysis of aircraft employ 
this method of counting. A crossing of a level with positive 
(or negative) slope is not made until the load has also 
crossed a second lower (or higher) preset level in the 


G2poo1be direction. Interpretation is hampered by 


a 





intermediate load cycles demonstrated in Figure 2. Two 
different load patterns are depicted but will produce equal 
counting results. However the loss of intermediate values 
is not a serious setback according to van Dijk [Ref. 5], who 
Senmcludema unat the majoricy of the load fluctuations of 
PMiLerestawere The larger excursions that originate from the 


steady-flight level. 


B. ACCELEROMETER LIMITATIONS 

The relationship between stress and acceleration at the 
Conver Obmelraviny 15 nol aS explicit as the relatiionship 
between stress and the strain measured at the critical 
point. The stress-acceleration relationship varies with 
weight, speed, altitude, and flight configuration. However, 
Lambert [Ref. 6] does not consider the variation critical 
for high performance aircraft since it can be assumed that 
all maneuver loads which will cause significant fatigue 
damage occur in the operational exercise of the sortie. 
Gust and maneuver loads that occur during the transit phase 
are assumed to be negligible by comparison. The stress-- 
acceleration relationship is considerably simplified by a 


Nroper choice of an operational weight and flight parameters. 


1.8 





IV. FATIGUE MONITORING PROGRAM 


The trend to extended operational service lives of mili- 
tary aircraft emphasizes the need for accurate fatigue damage 
calculations. Since protection against fatigue does not end 
with production delivery, subsequent responsibility includes 
further recording of service loads and environments, moni- 
toring and reevaluation in circumstances of mission changes, 
isDeectiOn ama mainvenance, an@istructural alterations 
throughout the service life. Due to the importance of 
fatigue strength in extending the service life of military 
agers) Vwmemomu. om Navy Hhasedevelopedgan extensive fatigue 
Nemicor nee rooranawhich ~wresentily monitors the individual 
structural fatigue status of 2700 aircraft. Expansion of 
the program is planned to monitor structural fatigue life 
expended on all Navy aircraft with the exception of transport 
anG@mmotary Wine aircraft. 

The Naval Aircraft Structural Fatigue Life Program enables 
StEplelLuralwinterrity decisidns to be made realistically in 
terms of fatigue experienced by the structure rather than 
by continuing the previously nebulous procedure of using 
PiemaUuMNbDe ry of Pen ele service or the accumulation of flight 
howe as #ammecepended liteser!terion. The Naval Air Develop- 
Demet eCtGei 1s cLaskea with Monitoring the structural fatigue 
life expended on each Naval aircraft. As a result a flecet- 


wide countings accelerometer prorram, combined with pilot's 





reports relating flight usage data, laboratory fatigue test 
Gata, and cumulative fatigue damage theory, was developed 
LO provide a measure of the actual Service life expended for an 


mic Toles aircraLt, throughous its service™life; 


A. UTILIZATION 
Usage of the fatigue data accumulated has provided the 

Navy with the source of information necessary to effectively 
aaiae mon laree si Meet of airerait., Ine ata are reported 
Oiamcverlyephe:. 6] and waseo eior: 

Early recogniticn of changes in service usage trends. 

imdenvarleavion and monitoring of aircraft accumulating 
unusually severe or excessive load factor occurrences and 
Ei rerun lew! Ualeually high fatipme damare eecumulation 
raves. 

Determination of aircraft time of retirement or 
rotation into or out of severe service uSage. 

Extension of service life based on changes in 
Micron flight restricticns, replacement or modification 
of major structural components, review of service history, 
agompertormance Of fulligsScale tavigue testse 

Provision of information effecting both organizational 
and in-service maintenance problems. 

Menizationset Gata wa the development/modification 


of naval aircraft design and test specifications/requirements. 
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B. DAMAGE PREDICTION 

AomWasepreviously mentioned, the main function of the 
structural fatigue program is to monitor the damage accumu= 
ae women Individual aiircraitw Typical of the method used 
to arrive at the damage calculation is the method employed 
to calculate the damage on A-7B aircraft. The following 
discussion will specifically apply to A-7B aircraft but is 
considered representative of the other types of aircraft 
monitored under the program. 

lie hull Seale Testing 

The first step in developing a damage calculation 

on the A-7B was to determine fatigue critical points. These 
points were analyzed and S-N data appropriate to these points 
were developed. A block program based on estimated loading 
per 100 hours of life time was used to fatigue test the 


aAinroeaun Vlable I), 


iu Negelde | 
A-7 Block Loading 


Maximum Cas) Number of load applications 
oad §(e8s} per load block 
Biljoekel1—4 Block 5-10 
2 Ns Acne) a7 OO 
SS So) 9D0 
2) 0 Sie 650 650 
4.55 450 N50 
Seas. BLG 250 
5 05 1 36 136 
6.65 Ny Hy 
(tas. aR a5 
8.05 yi u 
CaS al es 
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Based on thesproposed block loading program an 
estimated fatigue iife was predicted at each critical point 
based on Miner's cumulative damage theory. The full-scale 
fatigue test was then conducted (using blocks 1-10, which 
represent 1000 hours of flight time) until a crack initieted 
ateone of the critical points. The stress concentration 
factor of the S-N data was then adjusted to predict a fatigue 
iafe of 100 percent for the number of cycles that occurred 
On the full-scale fatigue test. The loads obtained from 
data reduction of accelerometer readings are then applied 
to the S-N data and a cumulative damage prediction is calcu- 
faved. Duel Gerthe uncertainty of the actual flight loading 
DuSGOrVeamsalelLymetactor Of two is applied to the damage 
pre daction- 

Aaeeeoda Data Reduction 

fie Benenerone ters used on the aircraft record the 
total exceedances at only four preset levels. The custodian 
squadrons report quarterly to NADC the accelerometer readings 
and pertinent flight data. Based on the four exceedance 
levels recorded, seven other acceleration levels are computed 
1Orevne  Losding spectrum fo be used for fatigue prediction 
in the case of A-7 aircraft. The computed exceedances are 
determined by use of the following equation [Ref. 8]: 


Cumulative Counts = Ae 


a 





Table II provides the values of the constants to be used 


for the seven levels. 


ALPE ey ir che 2 


A-7 Load Level Approximation 


Cumulative Counts A B @ 
5BDEc's il (Oly In(5-g counts)+tin(6-g counts) 2 
SE 2 ue IL G3 in(6-g counts)t+tin(7-g counts) 2 
7.5 gets 100 In€7=e counts)+1in(8~-g¢ counts) 2 
oe 1.00 21n(8=¢ counts)-1n(7.5-g counts) 1 
5a os .82 31n(5-g counts)-1n(6-g counts) 2 
0 ee 697 2in(4.5-g counts)-1ln(5=g counts) 1 
3.5 e's 95 e2ln(4-g counts)-1In(4.5-g counts) 1. 


The basis for these equations is obtained from reduced 
oscillograph data. An oscillograph recorder system is 
Insc Geom we number of fleevWaircralit which records time 
histories of airspeed, altitude, and normal acceleration. 
The recorder is energized by a switch on the landing gear 
retraction mechanism and consequently records data only when 
the landing gear is in the up position (the same is also 
true of accelerometer data so no landing information is pro- 
vided). Periodically these oscillograph data are compiled 
on aircraft [Ref. 8] and used to develop relationships between 
the various levels. The data are also used to redefine the 
stress-acceleration relationship since altitude and airspeed 
data are correlated with the loads. Unfortunately a semi- 


automatic data reduction system is employed which outputs 





data on all A-7's, and consequently, information on indivi- 
Gual aircraft is not available for a comparative analysis of 
trends between individual aircraft. The 11 levels of accel- 
eration now available are transformed into stress levels 

and the S-N data are referred to for damage prediction. The 
entire process has been adapted to a computer system which 
emavleis rapiid netrieval of pertinent Hlight and fatigue data 


Of any aircraft currently under this program. 


C. IMPROVEMENTS 

The fatigue monitoring program conducted by NADC has 
Mnauesticonably made “Significant advancements towards an 
accurate prediction of fatigue damage. There are possibili- 
ties for improvement based upon current knowledge of the 
fatigue phenomencn. 

It has been determined that maneuver loads occur randomly 
enawumeauep lock testing will not producesthe same fatigue 
damage as random load testing. Schijve [Ref. 1] mentions 
that, in general, the life in "equivalent" program tests 
fomuemecr Lhan in randomstests. ~Hesalludes!to ome Series 
of programs, where fatigue lives were about six times longer 
Lneoieiieirancomnm Load tesuse 

Through block loading fatigue testing on an A-7 aircraft, 
Miner's Rule has been altered by shifting the S-N data to 
account, somewhat, for the variation of load amplitude 
Peeourneut tite flight profiie, WhereiMiner's Rule inay be 
toe conservative, block loading on the other hand may be 


eve) (0) 0 malik Spa bow 
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The sequence effects of random loading cannot be 
ignored if accurate predictions are to be obtained. The 
remainder of this thesis considers the feasibility of 
applying sequential effects for fatigue analysis based 
on four-level accelerometer readings and the suitability 
fem pDivifigeeils Meunod Go a Computer. The next seerpaon 
will address the problem of developing a reasonable flight 
profile suitable for application to sequence effects and 


based on four levels of acceleration provided. 





ore DEVE COPMENT Or er fil Gin PROFLiE 


In order to apply sequential effects a reasonable flight 
Prouwele@inus L pe Produced thagesatisfiies Themiol towing 
requirements: 

1) The order in which the loads are encountered must 

be randomly distributed. 

2) The loading spectrum must be generated from the 

four levels provided by the accelerometer readings. 

3) Recognizing the individuality of aircraft flight 

profiles, two equivalent sets of accelerometer 
readings Must not Oe generate the same 
loading spectrum. 

Reference [8] and other similar reports from NADC con- 
Gaming Oslelllograph information on aircraft were analyzed 
iniecevetoOpins a flight profile’ suitable for high performance 


MiLiltGpamyeaircraft. 


A. OSCILLOGRAPH/ACCELEROMETER DATA 

Schijve {[Ref. 1] in Figure 3 differentiated between the 
gust spectrum and the maneuver spectrum experienced by air- 
Crom—elLne Cuctelogad SpeCcCLCrum COnSiISts Of angequal number 
of positive and negative small amplitude peaks occuring about 
Caematcadverlipht condition. ~The maneuver load spectrum 
Seem ausmeimthe remaining loads fan alrcratt experiences ana 
is highly asymmetrical about the steady flight condition. 


Lambert (Ref. 6] made a similar analogy as shown in Figure 4 
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and attributed the negative .loads in the maneuver spectrum 
as being due to "overswing" resulting from recovery from a 
high-g maneuver. From an extensive study of Swiss VENOM 
aircraft Branger [Ref. 9] concluded that the number of 
maximum peaks above the steady flight condition was equal 
to the number of minimum peaks below the steady flight 
Pendivion buL that the distribution is highly asymmetric. 
Oscillograph data [Ref. 8] support this phenomenon. A 
lower threshold of two g's in the positive loading regime 
and an upper threshold level of .5 g's in the negative 

min civeresiMempmoaqUuced ratiOs Ol positive to negative g's 
as high eas eight to one. The reason for this is demonstrated 
in Figure 3 where the number of positive loadings above 
threshola level A clearly exceed the number of négative 
loadings below threshold level B. Schijve [Ref. 1] indicates 
that all significant fatigue damaging loads originate from 
around the one-g level. Based on this evidence all 
accelerometer readings are assumed to be maximum peaks and 
likewise all oscillograph peaks in Bre pesitive Loading 
regime are also assumed to be maximum peaks Since the 
eriterion for recording @ peak is essentially the same 

as the accelerometer criterion. All loadings are assumed 
to originate from the one-g level or less (minimum peaks) 
depending, on the number of negative regime peaks the 


Oscillograph data provides. 
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B. ANALYSIS OF DATA 

Ditierent assumptions are made relating the peaks to a2 
common trend and the data are analyzed to verify the 
validity of each assumption. 

iP wml lativemioads Between Types 

Since all the aircraft considered were high perform- 

ance’ alircraitl, having missions that required high maneuvering 
loads, it was assumed that all aircraft experienced essen- 
tially the same number of loads based on 1000 hours of flight 
time. Figure 5 indicates the ToSsarithm ef Lhe number of 
loads that exceeded a particular g level for A-4E, A-T7A, 
F-~4B, A-6A and F-8E aircraft per 1000 hours of flight time. 
The data would not support the assumption that the maneuver 
loads of each type of aircraft were essentially the same and 
whus’ Maee mission requirements® of each type of aircraivt are 
not comparable enough to produce an equal number of similar 
l@adsamemowever, due to the similarity of the curves the 
relative ratios of load levels between different types of 
aircraft may be proportional. 

2. Ratio of Loads Between Types 

Although all the aircraft considered had missions 

agemanding high maneuver loads, some aircraft had missions 
Piao demanded mores loads#per flight than those of another 
Miecwatt ype. [In SspiteBol this te was assumed the loading 
spectra were similar and the ratios of loading levels 
between the types of aircraft are essentially equal. All 


. 


loading levels were ratioed to the lowest positive load 





level recorded. Figure 6 shows the ratio of five load 
levels for A-4E, A-5A, A-7B, F-4B, and F-8E to their respec- 
jiveme.0--e.59% oad counts. wo of the aircraft! s® ratios 
deviated considerably from the other three aircraft's ratios. 
Intermediate load levels were then investigated to avoid the 
possible effect of extreme gust loads at the lower level. 
Figure 7 shows the ratio of four load levels to the 4-4.5 g 
load level. The ratios again indicated a considerable degree 
of scatter. All loads that exceeded an upper level were 
ratioed to the load counts at lower levels. Figure 8 shows 
mic  rabwower suoullr 1Oads levels ve tne number of loads Una 
exceeded the four g load level. These ratios indicated even 
more scatter than those obtained in Figure 7. 

Since the load ratios investigateca did not indicate 
any clear similarity, it was concluded that the flight 
profiles between types of aircraft were distinctly different. 

3. Extreme Values 

Extreme values were investigated to determine whether 
there was a relationship between the load exceeded by extreme 
positive and negative loads. It was assumed that entry into 
these extreme loads was not intentional and probably flignts 
that encountered these loads encountered extreme values in 
both directions. Extreme values were defined as loads that 
accounted for one percent of the total number of loads 
occurring in either the positive or negative regime. Table 
J1I lists the number of positive and negative levels exceeded 


Gnewperceny ot the time. 
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TABLE III 


Extreme Value Occurrences and Level of Exceedance 


A/C Type Levey Now Eeceed: LE geal NiGwe Execedn 
A-HE 6.0 Beye Ae ~.41 SL ype 
A-6A Die 88.9 =.al 4.9 
A-7A Die Wit 3 -.1 20S 
F-4B oie: 42.1 ~ 41 Saye 
F-8E 6.0 LL - 41 35.4 


Analysis of Table III indicates no‘definite relation- 
Ship between the level exceeded in the positive regime and 
the level exceeded in the negative regime. A sample corre- 
lation coefficient of .642 for levels of extreme-value 
exceedance was not considered conclusive enough to support 
the assumption that a relationship existed between the 
positive and negative levels. Analysis of extreme values 
on a flight-by-flight basis is needed to confirm the 
assumption. 

ReLterence 10 conclaimam@se? | loprapneda Lo mol em oases 
"Biue Angels" Flight Demonstration Team compiled on a flight- 
by-flight basis. These data were also analyzed for a possible 
trend. By analyzing the data on an individual flight basis 
it could be determined whether or not extreme high loads 
tend to induce corresponding extreme low loads as a result 
of over-correcting recovery from the high load maneuvers. 
fable tv lists the nighést ang lowest = level recorded for 


feet laghtcs analyzed at random. 
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TABLE IV 


Blue Angel's Flight-by-Flight Extreme Values 


In IDs fed ole Maxie Ve 


1 ol Z 
2 pees ~ 
3 5.83 = 
4 16 = 
5 5 oloe . 
6 bat3 - 
ff Sue = 
8 5) 9 bss 
9 1230 - 
10 6.63 
i | 5.22 = 
12 6.49 = 


Detailed analysis of Reference 10 indicated that 
loads neither preceded nor followed an extremely 
and as Table IV indicates there is absolutely no 


[6oreassuming @ high probability of both high and 


gee 
Iss 
1.41 
Los 
a teyi 
slit 
105 
04 
mR es 
.08 
as 
1m Oa 


extreme low 


high 


load 


foundation 


low values 


Segurtredsom vhnegsame TlighuweOf particular significanecesin 


Warilesly is the lfact®that the lowest negative reading (=—2.67) 


Sceurred during a flight which hac the third lowest positive 


reading. It must be concluded that the negative readings 


occur randomly and independently of any extreme high loading 


ConGileLons »ssconsequentily, they could be considered as 


part 


of a "basic" Joading spectrum unique to the flying qualities 


Peermew articular aircraft consiaerved. The @xsteeme 


Gade 


low f¢ 





loads are probably due to high gusts loads and therefore 
the gust loads should be considered integrated into the 
"basic" loading spectrum. The dynamic "overswing" Lambert 
[Ref. 6] discussed must be above the .5 g level threshold. 
From the data investigated so far it appears that each air- 
Seaioehiaoma flieioy Loading protiile distinet from another 
type of aircraft and this profile includes all the negative 
g@ loadings and excludes all high g loads. If types of 
aircraft do indeed have a loading spectrum "fingerprint", 
ican Olly eapemoucyermined by comparing oscillograph data of 
individual aircraft of the Same type. : 
4, Individual Aircraft by Type 

In order to substantiate whether mee. of the 
Same type have identical flight profiles below some extreme 
positive load level, individual aircraft data are needed. 
A study conducted for NADC by Technology fncorporated [Ref. 
iiimecomcains individual aireraft oscillograph data on four 
A-6A's and six A+-7B's collected over a one-year period and 
representing data on 306.9 flight hours for the A-6A's 
anda. Jett sagemourscs tor the A-/Bes@. Although the number 
of aircraft of one type analyzed was small it was considered 
Sipmetmeals enouphevo decvermine af aircraft of the Same 
type do experience essentially the same loadin profile 
below some high load level. 

MEevaolvetcy sol Data 

The data were first compared with the oscillopraph 


Gata of Ref. 8 to assure comparable results between tne two 





reports. The A-6A data agreed quite closely and the A-7B 
data also agreed favorably in the lower loading regime but 
diverged significantly abcve the five-g level. This was 
NeoewcOnmsn dered Hmportant in this investigation, since the 
reeplonm of interest for developing a loading profile is 
contained below the five-g level. Since the ratios of loads 
wememcenstdered LO be vhne key vo developing’ a flight profile, 
the data were considered to agree favorably with other 
Scemme|@Goeraph davawsol the’same type aircraft. 
b. Comparisons of Each Type 

Mitac Mmmome ,ceedances fom anol vyidual alrerare 
of both types at each level were compared. The exceedance 
Esti nuloOmmror lndiviiaual gireratrt Gf each "type indicated 
Pen ume e rece Omecimllarity which would be necessary to 
develop a flight profile for each type. Figure 9 shows the 
number of loads exceeding an indicated level per 1000 hours 
Orgiliehte time forveach of the A-/B aircraft. The negative 
loads also exhibited a comparable Similarity. Analogous 
results were also obtained for the four A-6A aircraft. 

ite Muncie tsetse Occult Ine in each range were 
PoMOcCantGmeaneatoltrary level that mad @ range of ome pf. 
The one-g range level was chosen to agree with the ranges 
between recording levels of accelerometers installed in the 
aircrait. The number of loads that occurred within a prescribed 
range were determined by subtracting the number of exceedances 
at the higher level of the range from the number of exceedances 


at the lower level of the range. The-one gp range level wus 


Ur 
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increased and comparisons made until an unacceptable dis- 
agreement. resulted between individual aircraft. Acceptable 
agreement for the A-6A occurred only for readings up to the 
Bee Catiec we HOWever acceptable agreement for the A-7B 
occurred up to the 4-5 g range. Figure 10 shows the ratio 
of load occurrences to the number of 4-5 g occurrences for 
the A-7B aircraft. The load values were represented by the 
mean value of each range. The fact that one aircraft exper- 
lenced an increase in the ratio from the 2.5-2.9 range to 
the 2.9-3.3 range was not considered significant since that 
Bircralpmenadmerecorcea only 19e7 hours of flight time with 
the oscillograph installed. It was considered significant 
idiot se iC six dilirecrait had almost adentical values. 
e. Mh geht Profile 

Based On the above results ii was decided to 
eevertopea 1tlight profile for A-/B aireraft that would depend 
solely on the accelerometer counts in the 4-5 range. Table V 
indicates the load range and the ratio of that range to the 
4-5 g's each aircraft experienced. The ratio of loads in 
each range were assumed to be normally distributed about the 
ich lOm—DaseCmenmconeSCeaaba a Medan valve and a sample 
standard deviation were calculated. It was decided that a 
range of 1.65 times the sample standard deviation from each 
ratio mean would encompass an adequate range of loads that 


the aircraft would experience in each load level. Multiplying 
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A-7B Individual Ratios of Positive Loads to 4-5 g Loads 


TABLE V 


Load hoe ye SO 140.5 OM UIA Seal 76.9 
Range nes. . ies < Ines Ars. hrs. his: 
2.0-2.5 3.74 etal Beco 5.24 5 oes 2.85 
Pe o=—2.9 et: amen ie oo Baio fo eS! 
2.9-3.3 1.99 1 ee 46 ee. iL ee | 1.50 
3.3-3.6 1.14 Oe roag 1.49 ~714 855 
BO—3 29 13 .504 Raia ~ {70 5 ip) ay Os 
3.9-4.2 ANY Se 438 A414 400 439 
Wy2-4 5 ' ,298 Sos 5 SS ~ 304 2347 ~325 
4 5-4.8 Ser feos Mee e215 me | soo 
4.8-5.0 is) LAY £43 del pales Pe 


each ratio by the number of 4-5 ¢ readings would produce the 
actual number of loads the aircraft experienced for that 
flight period. 3.3 times the sample standard deviation pro- 
vides a probability of .95 that the aircraft experienced a 
load level in that range based on a Sample size of six. 
TaoVe Vi gives the range of ratios for each load level that 
the aircraft is allowed to experience in the positive regime 
for each 4-5 g¢ reading. 

mobples Vilveives agsimilar list of ratios: as Table 
ae L@y fie negative regime. The excessive number of negative 
readings aircraft (Serial No. 154370) experienced were not 
considered in the calculation of mean and sample standard 
deviation of the negative regime ratios. This aircraft had 


Only 19.7 nours of recorded joad time and its ratios were 
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TABLE VI 








A=7B Positive Ratio Range For Each Load Level 


Range Mean Low Mean High 
2s 2.399 4.01 55 OP: 
ee S06 1.92 EL Sys: 

: Pola es) 2 ie 
3.45 429 947 ae 
SIs 465 OSL .798 
4,05 538 rales 464 
4.35 »292 324 » 356 
4.65 mes .263 Pe oll 
4.9 lee L505 ~179 

JUN ey ici EIL 


A-7B Individual Ratios of Negative Load Range to 4-5 g Load 


Load Lom/ BO 24 140.5 VOI 113 4 faswrec 
Range hrs. hrs. Mesa lauatei esd Hise 
.6-.4 Iein 2) 7 9.02 4,807 9.24 6.74 By 8) i 
Hm 2 4.66 2.96 iL 2 Bjate 1.88 eee a2 .606 
-c-.0 ee 2 Oat . 3.6 408 Ses OSI 
ee oc) .274 oe 1a MG jOMAL .020 


not considered representative of the typical A-7B. Neglecting 
the excessive number of negative loads on this aircraft would 
not be more damaging in a cumulative damage calculation since 


. Saitple calculations with representative S-N data indicate a 
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loading with a larger range but a lower mean was less 
damaging than a loading with the same maximum peak that 
originated at one g and had a smaller range. Table VIII 
BiGeea Lec me lcm S BOL Fratton an the negative regime that 
have a probability of .95 of occurring for each 4-5 g 


reading based on a sample size of 6. 


ABE VIEL 


A-7B Negative Ratio Range For Each Load Level 


Range Value Low Mean Islabredel 
BSN. 5-91 (eS 5 10.34 
» 346 .260 1.673 3.09 
148 050 SSD, 667 
- .054 004 OS sy2(0:10 


Figures 11 and le display the range limits and 
mean for each load level for both the positive and negative 
repime,respectively. The mean value of each positive load 
level is considered to be an adequate representation of the 
foads"in that ‘range since the greater number of loads that 
occur below the mean value are offset in the damage calcula- 
tion by the fewer but more damaging high values in each 
range. 

5. sequence of Loads 
po comenmernerpreVvLOUssecOnsuraints a computer program 
was developed to arbitrarily assign a ratio for each load 


level to the 4-5 @ level reading. The value of each ratio 


55 





was determined by means of a normally distributed random 
number generator. The number of actual loads in each load 
level wiere then determined by multiplying each ratio by 
the number of 4-5 g readings that were inputted. 

lhnree Otherelevels of loads from A omsilemne renee? 
readings were also inputted that represented loads that 
occurred in the 5-6, 6-7, and 7-8 load ranges respectively. 
These load ranges were arbitrarily assigned mean values of 
5.75, 6.6, and 7.8 which were considered high enough to 
represent a conservative estimate for a damage calculation 
of a large number of loads occurring in this range since 
75 percent of loads recorded occurred below these values 
in each range [Ref. 11]. No attempt was made to distribute 
these high level loads throughout the entire range of their 
possible values due to the limited number of occurrences 
of these loads in Ref. ll. However if distribution functions 
for the extreme value loads were known it would be appropriate 
to subdivide the loads occurring in each range into a ame ee 
number of load ranges. Distributions similar to the one 
mused=in Taper il could be used but in this particular case 
phous c Eel puGlonm could noe be®subsStantiated from the data 
gee Lc emule. 

Based on the four load level accelerometer readings, 
all loads in the positive regime were assigned to be maximum 
peaks. Associated with the maximum peaks was an equal number 


of minimum peaks which were assigned values from the negative 








regime and the remaining number of required minima were 
assigned a value of one g. 

A uniformly distributed random number generator 
placed each group of peaks in a random order. Starting 
with a minimum peak, a load value was picked from each group 
ONanecl I feigllgl the Pe oem ee ales OUP miad been) Chosen...) A find i 
half cycle was then assigned to finish at a one-g level. 
Figures 13-16 indicate four flight profiles generated for 
accelerometer readings of three, five, three, and four 
representing loads occurring in the 4-5, 5-6, 6-7, and 7-8 
@ ranges respectively. 

It is highly unlikely that the order of loads gener- 
ated from the accelerometer readings would represent the 
actuai sequence the aircrart experienced during that period 
and thus a damage calculation based on this loading profile 
would be inaccurate. However, as the aircraft experiences 
an increased number of loadings over its fatigue life it is 
highly probable that the aircraft experiences a majority 
of the sequential loadings generated or ones quite similar 
to them since the loads generated are based on actual data 
obtained from aircraft of the same type whose order of 
ieacktinew 1S sown we occur randomly. The flight profiles 
generated will be of increasing significance during the 
HatLerewavares son the aircrain’s lite where the accuracy 
of a damage prediction calculation becomes of increasingly 


higher importance. 





Peon icaDl ling or Flight Profile 


The requirements stated previously for generation 
of a profile have been met. The loads are randomly distri- 
buted and are generated: solely from accelerometer readings. 
Two equivalent accelerometer readings do not necessarily 
generate the same loading profile. Due to the small number 
of aircraft sampled it cannot be assumed that the above pro- 
file is an accurate representation of the majority of A-7B's 
profiles. A larger sampling of aircraft would be necessary 
to further refine a representative profile. Currently the 
lowest accelerometer level being recorded on A-7B's is the 
five-g level. In order A adapt the flight profile developed, 
CoiemEevciovmle velerecorded would have to be lowered to the 
four-g level. The same is also true of the A-6A where the 
lowest level presently recorded is the four-g level and 2 
three-~-g level is considered necessary to obtain an indicator 
of the A-6A profile. However the above procedure demonstrates 
the feasibility of generating a flight profile from acceler- 
ometer readings since each type of aircraft apparently does 
have a "fingerprint" load spectrum below some high load 
level. The aiene prone Geveloped above will be considered 
representative enough of high performance fleet aircraft to 
proceed with accounting for the sequential effects ina 
fatigue damage calculation. The next section deals with 
the feasibility of applying a fatigue damage prediction that 
accounts for sequential effects and is suitable for computer 


analysis, 
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VI. ESDU CUMULATIVE DAMAGE HYPOTHESIS 


A cumulative damage hypothesis is needed that accounts 
hor ob Noor dc iemamn wate leads arewappliad to an aircrait that 
is suitable for computer analysis. 

The ESDU Cumulative Damage Hypothesis [Ref. 2] meets 
these requirements. Although originally intended for block 
loadings the method has been modified for analysis of loading 
on a cycle by cycle basis. The ESDU method requires no 
additional data over that required for use in Miner's Rule. 
Testing [Ref. 2] of lab specimens under block program loading 
has produced results that demonstrated a far better predic- 
tion of fatigue failure by the ESDU method than by Miner's 
Mile were nunber Of Similar specimens were tested until failure 
which on the average took 30.9 programs to occur. : The Miner's 
Rule estimate for failure was ll./? programs; whereas the ESDU 
estimate was 28.7 programs, an improvement of over 240 
melee nin. 

The basic hypothesis of the ESDU method is centered in 
the concept that. once the stress applied to a component ex- 
Cece CmpmolOtmy mona) limi, restdual stressesgbuild up 
Which Be reac Maat er the actual mean stress of subsequent 
stress cycles. The following basic assumptions are made 
when applying the ESDU method: 

1) The material is perfectly elastic below the 
yield stress and ree plastic above it. The yield 
Sete osmeee Lower as tiemmeem percent Offset stresis. 
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2) The stress-strain characteristics of the material 
are the same in tension and compression. 

3) The stress concentration is constant for all 
Momiiawlmotresses . 

4) Yielding is localized to stress concentrations 
and therefore does not significantly affect stresses in the 
eocoeeOlLecne member. 

5) A constant amplitude S-N curve for the material 
is available. 

A member that is Subjected to a nominal stress whose 


peak value is SMAX has an elastic stress at the stress 


concentration of SMAXEL equal to: 


SMAXEL = KT +» (SMAX) (2) 


If the value of SMAXEL exceeds FP, the .2 percent offset 
yield stress, the mean value of the subsequent cycles will 


be lowered by an amount equal to H where: 
H = FP ~ SMAXEL (4) 


mecCeocivicwmmio Aloo lruc wor am GxXCeESSive compressive stress, 
which QuemeeeLnNeCeDUald wp Of residual tensile stress causes 
an effective raising of the mean stresses of the subsequent 
Pyelcomemtidenres 7-19 indicate the resulting effect that 
the different types of loadings have on subsequent cycles. 


In Figure 17 neither the maximum nor the minimum peak exceed 
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FP so no alteration of the mean stress in subsequent cycles 
is required. Figure 18 indicates the effective raising of 
the mean stress of subsequent cycles due to a compressive 
stress that exceeded -FP. Figure 19 shows the effect of 
compressive residual stresses due to a tensile loading | 
Ricceomia Lm value exceedeq ThE | 

Application of the ESDU method is quite simple and is 
easily adapted to a cycle-by-cycle analysis. <A possibie 
objection to the validity may be that the effects of stress 
relaxation are neglected, which causes the residual stress 
to fade away if given enough time. However Schijve [Ref. 1] 
states that this effect has not been observed in aluminum 
alloys, the materials most aircraft wing components are 
made from. 

The next section discusses the computer program used to 
develop a flight profile from accelerometer readings and to 
apply the ESDU Cumulative Damage Hypothesis for a damage 


Pemmoulation of this profile. 
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~ VIL. DISCUSSION OF COMPUTER PROGRAM 


Program ESDU was developed primarily to demonstrate the 
feasibility of adapting the ESDU Cumulative Damage Hypothesis 
VOma computer program for use’ in calculating the fatigue 
damage of high performance aircraft with accelerometer 
readings as data input. In order to provide a concrete 
example the program was specifically designed to provide 
damage calculations on A-7B aircraft that had accelerometer 


recording levels at the four, five, six, and seven g level. 


A. FEATURES 

There is virtually no limitation to the number of readings 
at each level (ten digit maximum); however, only a limited 
number of readings will be used in the calculation if the 
number of four-g level readings exceeds Seon NCE CrOogran 
requires 152 K bytes of storage and less than three seconds 
of computer time for each aircraft calculation. Data input 
consists of the number of aircraft to be analyzed, the 
serial number and four accelerometer readings for each 


perrcralt. 


B. PROCEDURES FOR CALCULATION 

The program is essentially subroutine structured with 
the main program linking the ten subroutines required for 
imaeecaleulation, ‘The’ main program accepts the input data, 


assures compatible storage requirements, requests maximum 
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and minimum peaks, requests a random ordering of the peaks, 
Gmetad hy ea i Seormamcalculation of damage for the light 
sequence using the ESDU method. A brief description of 
Som suotrourine L£Ollows . 
i our OUtarte” Fix 

Fix assures the storage requirements will not exceed 
152 K bytes. If moré than 325 four-g level readings are 
PHimeneseOdm=Lor One alrerait, 1t is hypothetically possible to 
generate more than the 10,000 peak storage limitation. If 
more than 325 four-g level loads are inputted the four 
recorded levels are ratioed so that 325 four-g level loads 
Will be provided to the main program. Fix also provides a 
correction factor to compensate for the reduced number of 
peaks in the final damage calculation. 

Zee oUbcOouULne Upspec 

Upspec provides the main program with the number of 
maximum peaks and their values to be used in the damage 
calculation. A ratio of each load level below four g's to 
the number of 4—~5 g loads is determined by a normally 
distributed random number generator provided by subroutine 
Gauss. If the random number provided exceeds the .95 
percentile range for each load level a new number is 
requested. After the number of loads in each level are 
determined, subroutine Hiload assigns the mean g load value 


to each peak in each load range. 
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3. Subroutine Dnspec 2 


Dnspec performs essentially the same function as 
Upspec, except it provides the main program with minimum 
peaks. Appropriate ratios are assigned to four levels of 
minimum peaks. NPCoamee erination of the number of peaks 
required in each load level the remaining required minimum 
peaks are assigned to a fifth level, which represents the 
one g minima, to equalize the number of maxima and minima. 
subroutine Dnload assigns the load value to each of the 
five ranges. 

4, Subroutine Mixer 

Mixer provides the main program with the maxima and 
iiniemarranced dned random order sO as to represent a flight 
profile. Mixer requests Subroutine Shufle to randomly 
arrange both the maxima and minima py means of a uniformly 
distributed random number generator provided by subroutine 
Random, Mixer then alternately picks a value from both 
groups and stores them in an array mae] all peaks have been 
used. 

Sew mOlt Lie Molu 

ESDU calculates the damage incurred vy the aircraft 
on a cycle-by-cytle basis using the ESDU Cumulative Damage 
Hypothesis. ESDU requires data about the particular aircraft 
being analyzed. <A stress concentration factor, S-N data, 
ultimate tensile stress, .2 percent offset yield stress, 
Mma@mecOnversl Om TrOMmss tos Suresseare neededifor thel aircraft 


Pemeie ae eie wer pi1cal Point. Folynomial approximations 
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We vPemiccim naw DOD gata bub more™precisile descriptions 
of the S-N data would be desirable. Linear interpolation 
of tabulated S-N data would be an acceptable method. S-N 
data for only one mean value are all that is necessary 
monde *the program AGG UBOS ie = data tc accommodate 
different means by ane method described in Ref. Zee) onc 


next section sample data runs are discussed and the 


relevancy of their-results to actual aircraft predictions. 
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VIII. DATA ANALYSIS 


Actual accelerometer readings reported from fleet air- 
craft during heavy maneuver periods were fed into the program 
to determine if reasonable damage predictions would result. 
Although it was realized that the damage predictions have 
no realistic value except for A-7B's due to the flight pro- 
file variations between types of aircraft, subroutine ESDU 
was altered to accommodate the pertinent load data for A-6A 
and F-4B aircraft as well as A-7B aircraft. A fourth "hypo- 
thetical" aircraft was devised in order to cover a range of 
Sicess COMNCeNLatlon Lactors and @-tc-stre'ss conversions 
eommonly found in high performance aircraft. The stress 
concentration factors and g--to-stress conversions represent 
the actual values for the fatigue critical points of the 
aircraft types tested. The g-to-stress conversion is deter- 
mineca by the design limit load of the aircraft, which is 
the amount of nominal stress the fatigue critical component 
experiences at the 6.5 or 7.0 g load depending upon the 
aircraft in question. The design limit loads varied from 
HOvOOmpe simon the A-6A to 24,600 p.s.i. on the?A-7B. 
Pelecata was obtained DY interpolation of *Goodman Diagrams. 
contained in Ref. 12. The material used at the fatigue 
Porm@ioalempointe for All eGireraft was ‘aluminum plate 7075-T6. 
Table IX provides the pertinent information on each aircraft 


used in the test runs. 
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TABLE IX 


Aircraft Tested 


Aircraft F-4B A-7B A-6A Hypo. 

Ke Syees: Se 3.0 ange) OS 

stress (KST) 

per g 4.723 3.5143 5) 5 dU! 57 

Location outer outer inner 

eriitwmeal@et. wing wing wing = 
panel panel panel 

Flight hours 25.8 31.5 33.4 cules 


Five test runs were performed on each aircraft to deter- 
mine if the prediction variation was significant with differ- 
ent flight profiles. a for the 4-5 g levei of the A-7B 
were developed from the equations provided in Table IIl. 
mache LOCiLi1 Cale | Lent prom le was Subjected to ai Miner's 
Rule approximation for comparison with the ESDU method. 


Tables X=-XIII show the results of these runs for each aircraft. 
TABLE x 
P-4B Test Runs with 25.8 Hours of Combat Flight Time 


Meceleracticonsdata (2's) Herne oredr err on 


Four Puigve Six Ioana ESDU Miner's 
30 ial 2 0 0005624 POM COG 
30 ital 2 0 .0006568 Sealey s) S10 
30 aga Z 0 | .0008384 isa 
30 SMa 2 0 .0009822 5 OES y fee 
a0 11 2 0 0010484 -Bilsisreleytt 
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TABLE XI 


A-6A Test Runs with 33.4 Hours of Combat Flight Time 


Accelerometer data Damage prediction 


Pour Five Six Seven ESDU Miner's 
26 2. 0 0 OOO L st .0147439 
26 5 0 0 .0012068 5 (oaL 9) TeSys: 
26 5 0 0 CGH 563 OlGeeg 
26 5 0 0 .0012042 SONS 67/6 
26 5 0 0 .0014054 OMG ay 

TAB CI 


Four 
194 
194 
194 
194 
194 


Accelerometer data 


Five 
120 
120 
Let 
1210. 
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Sale 
Di 
on 


57 


5 
Bf 


Seven 
ie 
12 
12 


i 


Damage prediction 


ES DU 


.0018248 
0017034 
.0017585 
.0017833 


mooG 320 


n= (pebest hums with 28 ..5eHours of Combat Flight Time 


Miner's 


.0118987 
Enno c 
ROM fe3 
.0119610 


POMeG 3.6 





TABLE XIII 


Hypothetical Aircraft with A-7B Flight Data 


Acceleration data in g's Damage prediction 

Four Five Da seven ESDU Wik igesin 
194 120 Bul, lia .0016852 .0058465 
194 120 yf) 12 .0014576 ~ SHOE SSA, 7 
194 Zhe Sif ges .0014858 OMI Ey 
194 nae ai ie .0013641 nea oui fal 
194 120 ot 12 0012437 [0053356 


In order to easily interpret the degree of variation 
peEween ioicmeves tl lighe protates s the mean value of the 
damage calculation and the percent of deviation from the 
mean of the largest excursion for each aircraft was calcu- 


lated. The results are tabulated in Table XIV. 


TABLE ¥IV_ 


Results of Test Runs 


Aircraft P-4B A~-7B A-6A Hypo. 
ESDU mean .0008216 .0017404 .001196 .0014473 
Percent excursion 29.11 4.847 ass 16.44 
Miner's mean FO15596 Owaiy 3 OM Sif 9 0059423 
Percent excursion eee 55 ome BED 
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some definite conclusions were made from the test runs. 
The ESDU method makes a considerable difference in a damage 
prediction as compared to Miner's rule. The variation from 
the mean in the damage prediction was noticeably larger in 
the ESDU method than from Miner's calculation, although they 
both were based on the same flight profile. Apparently the 
Varilaglon im, tne mumber of Joads in the profile is not quite 
fom ose readnl sep nes damage caleultavion as is the order of 
the loading since the ESDU method, which accounts for the 


MieletmoLe Gic Toads, Nad 4a much Parser deviation. 


Pee LILBE ES TaLMA TE 

It was decided to calculate the life expectancy of an 
A-7B based on the number of occurrences of loads per 1000 
hours from the oScillograph data of Ref. 8. Due to the 
large number of loads involved the data were subdivided 
into blocks of 26 representing 38.46 flight hours each. The 
loads were also divided into blocks of 52 representing 19.23 
flight hours each. Ten runs were made on each set of data 
and by using the mean damage value each damage calculation 
was multiplied by the respective number of blocks to produce 
the damage accrued for 1000 hours. The Same procedure was 
fomr1Oweo nor ae Miner*s Rulle prediction. The results are 
Shown in Table XV. | 

Increasing the number of blocks from 26 to 52 resuited 
in an increase in the damage life for 1000 hours of 17.9 
percent whereas the Miner's damage prediction showed a 


negligible increase of .396 percent. ‘The reason for this 
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TABLE XV 


A-7B Damage Accumulated for 1000 Flight Hours 
‘Based on Accelerometer Readings of 146, 
Gee Lom, Snow) ¢or use Bilocks 


Method No. of Blocks Damage/Block Total Damage 
ESDU oa 0008126 021128 
ESDU ae .000479 .024908 
Miner's 26 .00405351 . 1053913 
Miner's ve .0020348 .1058096 


imoereasc@its due to the factt that everytime agblock begins, 
the material is essentially free of residual stresses, and 
thus no favorable mean stress correction occurs until a high 
EB load causes plasitic deformation. The Miner's approximation 
is net concerned with the build-up of residual stresses, and 
thus should show no significant change in damage prediction 
as the number of blocks increases. 

It was assumed that after every flight the favorable 
residual stresses were relieved by the damaging effect of 
landings and the taxiing phase. It should be noted that the 
taxiing phase will only have a damaging effect on fatigue 
critical points outboard of the main mounts. In addition it 
was also assumed that the Gailesespredicusen increaseq Tinearly 
as the number of blocks increased. Based on these assumptions, 
a damage calculation for the 1000 flight hours was divided 
into 500 flights that would represent an average of 2 hours 
between actual residual stress relieving. The aamagre 


Calculated based on a.linear relationship resulted in a 





Gamaze premicpicneot .09 or an expectedi life of 11,000 hours. 
This was not considered totally satisfactory since NADC's 
predictions are producing life estimates [Ref. 4] in excess 


of 22,000 hours (excluding the factor of safety). 


Eee FECTS Of STRESS RELIEF 

lt was decided to investigate the AE Remains between 
increasing flight blocks and damage prediction. The A-7B 
previously used for a one-block flight was divided into 
flight blocks of one, two, three, four, six, and twelve. 
The damage prediction was then calculated for the entire 
fioenourss@tetiligcht time. Phe results, including Miner's 
prediction, are indicated in Table XVI. Ten runs with each 
set of data were run with the mean value of the damage for 


the ten runs tabulated. 


TABLE XVI 


A-7B Damage Predictions with Increasing Flight Blocks 


No.eet) Blocks Damage/Flight Total Damage 
ESDU Miner's ESDU 
il .OOLTHO4Y nO / (oe 0017404 
2 .00106697 gOS ee BCC chs sy 
BS} nO OO (O07 5 .0019376 pOO cme 2S 
4 .00058844 01185008 00235376 
6 .00043121 mode sles .00025873 
lee .0002321.84 POW? 6695 (Oe 7 ooe 
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Except for the decrease in prediction from a block of 
three flights the damage is clearly increasing as should be 
expected. The predictions from Miner's Rule are relatively 
Stable. Although the increase in damage prediction is not 
Pricarivereclated CO an increase in flight blocks, 
the relationship indicated a linear increase with the 
feeorrunm Of =the@number of ®flight blocks, ~Figure 17 indi- 
cates the trend. Based on these results an aircraft damage 
prediction could be determined by using a one flight block 
run and an additional block run that would be compatible 
Monies deve. hen assuming a linear relationship between 
imiemoamertc amd the logarithm ofgflight block, the damage for 
any number of flights during that block could be determined. 
Based on this relationship the damage home tne, 500) minehec 
plockeOm the Csicillogpraph data ior 1000 flight hours produced 
a damage prediction of .0372 per 1000 hours or a life expec- 
tancy of 26,851 flight hours. This may seem a rather high 
estimate,but the loading profile the aircraft was designed 
for is considerably higher than the actual loads experienced. 

igo tsemey supfestedsthat the Jinear relationship of the 
Boca ienn sOnerhe Mumnber of PTlirnts may clearly a linear 
relationship over the entire range. Extending predictions 
based on 26 or 52 flights for 500 flight blocks would not 
produce accurate results. However if the number of flights 
was close to the blocks used in obtaining two point predic- 
tions the linear relationship should produce reasonable 


Pecuwlgemeum tne ease Of the A-/B with 31.5 hours of 


com) 
" «4 





accelerometer data,blocks of one and twelve flights should 
predict a reasonable damage prediction for the aircraft if 

Mo ace elem a 1lfhntsmwere involved sduring this) reportine period 
which results in an ae eee 1 sng Lc perm lI ShiGwOm 2 ai 
mOours. 

How much relief of residual stresses is caused by the 
landing and taxiing portion of the flight is presently 
unknown. Assuming that all the residual stress is relieved 
as a result of the ground phase is as unreasonable as assuming 
that none of the residual stress built up by a high load is 
relieved. Somewhere in between these two extremes is the 
fe cloi anOUmmuotescpress melliei that Occurs. Until this 
information is known the assumption that total stresis relief 
eccurs aiter each flight Ws a highly conservative estimate 
that would produce damage predictions that have a built-in 


factor of safety. 


56 





IX. SUMMARY AND RECOMMENDATIONS 


A need has been demonstrated for the implementation of 
a damage hypothesis in fatigue prediction that accounts for 
sequential effects. Predictions based on block loading 
programs have shown to be an optimistic prediction of the 
fatigue damage. Recording devices that have a direct rela- 
Prous ip sem bne SUress e€xperi1enced at a fatigue critical 
point such as time history strain gauges are highly desirable. 
Due to the nature of aircraft loadings, most fatigue damaging 
loads are initiated from the one-g level or slightly below 
it. As a result of this type of loading, accelerometers were 
considered acceptable for indicating the number of loads that 
occurred above the preset levels. Unfortunately no informa-~ 
tion was provided concerning the order in which the loads 
SecuUrreaw elt was found that the occurrence of a particular 
load in flight could not be sequentially associated with the 
occurrence of any other load. As a result it was concluded 
that maneuver loads occur randomly. 

Based on accelerometer data, a flight profile was devel- 
oped that includes all fatigue damaging loads that a 
Banptcicular type aireratt experiences. Unfortunately the 
load levels now being recorded on the A-6A and A-7B aircraft 
are above the level necessary to obtain the flight profile. 

The need for recording load levels above the eight-g 


level in A-78 aircraft is appreciated due to the highly 





damaging effect these loads impose; however, analysis of 
flight data indicates that loads above the eight-g level 
occur on the average less than five times for every 1000 
hours of flight time. Reference 13 eee that loads 
with a peak of 8.05 g's produce a damage of .0001. Based 
on these data an A-7B which experienced only 8.05 g's could 
iO eoree,0O,0008nours belore Lfatvigue failure. It is 
Beeconmendeadeuvnat Ghe accelerometer levels on A-7B aircraft 
Dew rowered GO anclude readings at the four-g level. 
Estimation of eight-g occurrences could be provided by 
allowing a reasonable number of the seven-g loads to occur 
above the eight-g level by means of some distribution 
Similar to the equations provided in Table"III. 

honGaiaeT Me wOraer Oust liging loegs has 4a signaficant 
effect on the damage calculation provided by the ESDU method. 
Misti =enCcman.LeC spar sen spiemaircralg, 16 18 anticipated that 
the majority of the loading sequences generated would have 
actually been experienced by the aircraft. Varying the 
oa VOmemelOaddSs In ach rangeshad 4 leissisignitacant eitfect 
on the damage prediction, since the Miner's Rule approximation 
produced approximately similar results regardless of the 
MumMDeEE Clem Oaus a2 particular trial run peoauced” Neglecting 
a variation of load level ratios would preclude the possibility 
of one aircraft accidentally being assigned an extreme low 
rit Vomuor Cache nlaghny period. 

The flight profile generated from a small sample of 


aircraft cannot be considered as conclusive evidence that 








each aircraft has a "fingerprint" flight profile below some 
high load level. Oscillograph data now compiled by type of 
Per ecreteeanOuld De requced Moimpan individual aircraft basis 
so that the profile hypothesis can be substantiated, 

The ESDU cumulative damage hypothesis provides a reason- 
able explanation of the overconservativism of the Miner's 
approximation. The ESDU method is suitable for computer 
analysis of damage prediction of a large fleet of aircraft 
due to the short run time required for calculation. 

Assuming that the residual stress imposed on the air- 
craft component due to a high load would remain indefinitely 
imUnccasOvamlcwmmnoweyer Lhe amount of stress relief caused 
by landing and taxiing loads is not available. It is recom- 
mended that the amount of stress relief imposed by the ground 


mease Of the flignt De’ siven considerable attention to enable 


Pe licuedse@maecUlracy or damage prediction by the ESDU method. 


Witil information is known concerning the ground phase of 
the flight, each flight could be assumed to relieve completely 
the favorable residual stresses that developed during the 
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FOURsFIVE,SIX, SEVEN, PREDIC 


NPEAK, PREDIC] 


JST 
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F4.40,7T29; 


9 


NACRE! *,%14,*FOUR® .7129,°FIVE® -1T445*SIX*¢ 
S) 


Veae PRED Clie 7/4 


STOP 
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M@O@UGerlrcCeuClUeCCOCCOCCCC CCCCCCCCCCEGCCCCCCCCCCCCCCCCECCCCCCCE 


CCCCCCCCCCCCCCCCCCE 


C 
C 
C 
C 
C 
C 
C 


SEQUES(10001) 


DNPEAK(500Q1; 


(5;¢)} NUM 


IMENS IGN UPPEAK( 5000), 
ITE (6,6) 


- 


D 
WR 
RE 


tO 
(FOUR,FIVE, SIXsSEVEN,AJUST) 


) NACRFT, FOUR,FIVEsSIX,SEVEN 
T2325.0) GO 


Ousm OC Ot 
OY TOMO 


— 


G PEAKS 


ASSURES MEMORY REQUIREMENTS SATISFIED 


UBROUTINE TO GENERATE MAXIMUM 


S 
A/C HAS EXPERIENCED 


OVQVOVOOO 


eNACRFT,FOURsSFIVEsSIX, 


GENERATE MINIMUM PEAKS 


" LL DNSPEC (NPEAK,NACRET, ONPEAK; FOUR; IIIT} 


WOU 


A 
F (NPEAKeEQ-0) GU TO 4 


INTEGER FQR 


D 
RATOK 
TFIXC FIVE +IFIX(CSIX)+ 


OOOO 


ORDER 


3 CALL MIXER (UPPEAK, DNPEAK, SEQUES,;NPEAK,NODD) 
EATIGUE DAMAGE 


NPEAKyPREDIC) 


LOADINGS PLACED IN RANDOM 
T 
’ 


OBTAIN 


OOO YOOYO 


FOURsFIVEsS1TXsSEVEN, PREDIC 


oh OP) nd 


O 


STOP 


momOUrwsicS;s*FIVE' , 1445 ° SIX"; 


{/) 


“4. 
f 
6,7139+F4.0,729;F4.03;7T%433F4.053;7T60,F4.0; 
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mnie (Ombre COOCOCCCCOCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


C 
C 
C SURRUUISENE MaEDUCES DAtA FOR COMPATIBLE 
: MEMORY STORAGE WITH DIMENSIONED ARRAYS 
& 


Vere wUMMeEmeeeeeeOeOCCeCeCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


SUBROUTINE FIX (FOUR,FIVE,SIX,;SEVENsAJUST) oe 
X=325.0/FOQUR 

FOUR=325.0 

FIVE=X*FIVE 

Sie Ges 10 

SEVEN=X*SEVEN 

IF (AMGDCFIVE,1.)-GT.0.5) FIVE=FIVE+#1.0 

NP CAMGO(SIXG1e).Gi.0.5) SIX=SIA+1.0 

IF (AMGDESEVENs Le} eGT.025) SEVEN=SEVEN+F1.0 
FIVE=AINT(FIVE) 

SIX=ALN@pS 1x) 

SEVEN=AINT(SEVEN) 

AS t= be or x 

RETURN 

END 


ea ee ea ce oe CCE CCE CECCECCCCECE CELE 
E SUBROUTINE GENERATES MAX G PEAKS 
PPC COmUrCCOCMmemeeCeCUGGeCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


SUBRCUTINE UPSPEC (NPEAK,NACRFTs; FOUR, FIVEsSIXsSEVEN; a 
LUPPEAK; IIIT) 

c OIMENSITON UPPEAK(5000)—, NPEK(12) 

© ae se NORMAL DISTRIBUTICN CHOOSE NUMBER 

OF PEAKS IN EACH LEVEL 

c 
DO 1 f=1;,12 

1 NPEKCI) =0 

: 

2 NEED OOD INTEGER TO INITIATE CALL 
NOGDD=1LFIXCFOURIFIFIX(CFIVEtIFIX(SIXI tT FIXCSEVEN: 
teen. EO.OI Camm Ue ls 
NCDD=(NACRFET/NODDI*IITI 
IF (MOD{NOCDO,2)-EQ.1) GU TG 2 
NOSDD=NGOO+1 

A Z2 ITA%=NOD0 | 

, FIRST RANGE 

Z Pome s( FOURPSEO.O) GO TO 12 

C NORMAL ODISTIKBUTION RANDOM NUMBER 

é GENERATOR 

3 CALL GAUSS (1%,.9754,4.91;V) 
PeeGvGitsreae) GU 10 4 
ittevee tree 399) GO) 10.3 
V=FOQUR*V 
NPEKCLISIFIXOV) 


OY 


POUNCE Ort STO UNERREST WhoLE NUMBER 
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SECOND RANGE 


IF (AMOD(Ve1l-.)-GT.005) NPEK(1) 


OO OOO 
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